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ABSTRACT 


The  numerical  f rontogenesis  model  of  Williams  [1974]  is 
converted  to  in   p  vertical  spacing  and  modified  to  include 
a  higher  upper  boundary,  tropopause  and  stratosphere,  mois- 
ture and  variable  vertical  turbulent  diffusion  coefficients 
of  momentum  and  heat.   The  numerical  solutions  show  more 
realistic  quasi-steady  fronts  forming  within  1-2  days. 
These  solutions  are  examined  and  compared  over  a  range  of 
various  boundary  value  coefficients  and  various  ranges  of 
temperature . 
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I .   INTRODUCTION 

Hoskins  and  Bretherton  [1972]  have  shown  analytically, 
and  Williams  [1967,  1972]  has  shown  numerically  that  dis- 
continuous fronts  can  form  within  a  finite  period  of  time 
if  no  turbulent  diffusion  is  present.   These  studies  sug- 
gest that  a  discontinuity  will  form  within  24  to  36  hours 
when  reasonable  initial  conditions  are  used.   If  turbulent 
diffusion  is  present,  it  can  be  expected  that  a  balance  will 
be  achieved  between  the  f rontogenetic  advections  and  turbu- 
lent diffusions  of  momentum,  heat,  and  moisture.   The  front 
should  remain  in  this  state  of  quasi -balance  as  long  as  the 
large  scale  deformation  field  causes  f rontogenetic  advec- 
tions around  the  front. 

The  numerical  f rontogenesis  model  of  Williams  [1974] 
shows  the  effects  of  horizontal  and  vertical  turbulent  dif- 
fusions of  momentum  and  heat  in  the  formation  of  these 
quasi-steady  fronts.   This  model  which  includes  an  Ekman 
boundary  layer  will  hereafter  be  referred  to  as  W74  and  the 
model  as  it  existed  in  Williams  [1972]  will  be  referred  to 
as  W72.   Cornelius  [1974]  added  moisture  with  its  subsequent 
condensation  and  release  of  latent  heat  to  W74.   This  in- 
vestigation raises  the  upper  boundary,  includes  a  tropopause 
and  stratosphere  and  appends  variable  vertical  turbulent 
momentum  and  heat  diffusion  coefficients  to  the  work  of 
Cornelius  [1974]. 
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The  purpose  of  this  study  is  to  obtain  and  examine  frontal 
solutions  throughout  a  four-day  period.   The  physical  model 
is  essentially  the  same  as  one  of  the  models  treated  by 
Hoskins  [1971]  and  Hoskins  and  Bretherton  [1972].   The 
basic  model,  as  it  now  exists,  numerically  describes  frcnto- 
genesis  which  is  forced  by  a  nondivergent  horizontal  wind 
field  which  contains  stretching  deformation.   This  deforma- 
tion wind  field  is  constant  in  time  and  independent  of  height, 
except  in  the  Ekman  boundary  layer.   The  moist  hydrostatic 
primitive  equations  with  diffusion  are  used  in  which  the 
time  dependent  quantities  are  functions  of  y  and  p  only. 
Therefore,  output  from  the  model  is  in  the  form  of  vertical 
cross -sections  in  the  Y-P  plane. 

The  model  of  W74  has  been  transformed  to  pressure  coordi- 
nates to  obtain  more  accurate  values  of  pressure  for  the 
moisture  equations  of  Cornelius  [1974] . 

In  Section  II,  the  basic  forecast  equations  are  developed 
and  a  simplification  which  keeps  the  problem  two-dimensional 
is  discussed.   The  initial  conditions  are  presented  in  Sec- 
tion III.   Convective  adjustment  processes  and  latent  heat 
release  in  the  model  are  included  in  Section  IV.   The  experi- 
ments undertaken  are  examined  in  Section  V  and  conclusions 
are  given  in  Section  VI. 


13 


II.   BASIC  EQUATIONS  AND  BOUNDARY  CONDITIONS 

The  hydrostatic  primitive  equations  with  diffusion  and 
moisture  may  be  written  as  follows: 

9  V  9  (o)V) 

_  +  v-(VV)  +  ___  +  v<f>  +  f(kxV)  =  AMV2V 

a     ^ 
+  W     (KMFp')  >  C2.1) 

||*  v.(ev)  ♦  y|  (.9)  -  Aev>e  .^  (Ke||) 

+  QA  +  QM  +  Hc  ,  (2.2) 

|£  +  V-(qV)  +  ~  (coq)  =  A  V2q  +  |-  (K  I3-  ) 
9t       ^—^    dr>         ^  Q  3p  v  q  9p  J 


*   MA 

"    MC    - 

v.v*|£-  o, 

d4  -  - L<w  > 

where 

L(p)    =  £  (  E-  )K  = 
P       P0 

9 
CP    3p 

(2.3) 
(2.4) 
(2.5) 


(  ?-)*. 


The  equations  are  in  pressure  coordinates  with  slight  modi- 
fication due  to  the  addition  of  the  equation  of  conservation 
of  water  vapor  (2.3)  which  is  analogous  to  the  thermodynamic 
equation  (2.2).   Sinks  and  sources  of  heat  and  moisture  have 
been  added  to  the  thermodynamic  equation  (2.2)  and  the  con- 
servation of  water  vapor  equation  (2.3). 
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The  quantities,  A.,,  A„  and  A   are  the  horizontal  turbu- 
lent diffusion  coefficients  of  momentum,  heat  and  moisture, 
respectively,  while  K„,  K   and  K   are  the  corresponding  ver- 
tical coefficients.   These  coefficients  were  assigned  various 
reasonable  values  in  the  following  experiments.   Dry  and 
moist  convective  adjustment  processes  are  represented  by 
the  functions  Q.  and  QM  while  latent  heating  due  to  large 
scale  condensation  is  represented  by  H  .   The  quantity  M. 
is  a  sink  or  source  of  moisture  due  to  moist  convective 
adjustment  in  the  moisture  equation,  while  NL,  is  a  sink  for 
moisture  due  to  large  scale  condensation. 

Boundary  conditions  for  the  above  equations  are  such 
that  the  domain  has  a  flat  lower  boundary  and  a  constant 
pressure  surface  where  go  =  0  as  the  upper  boundary.   The  no 
slip  condition  is  used  at  the  lower  boundary  and  the  no 
stress  condition  is  used  at  the  upper  boundary.   The  condi- 
tions are  as  follows: 

w(x,y,0,t)  =  w(x,y,pt,t)  =  0 
V(x,y,0,t)  =  0 


8V 

km  3p  Cx,y,Pt,t)  =  0 


Ke  |£  (x,y,0,t)  -  KQ  §i  (x,y,pt,t)  =  0 

Kq  y£  (*,y,o,t)  =  Kq  |a  (x,y,Pt,t)  =  o 


(2.6) 


wh 


ere  p   represents  the  pressure  level  at  the  top  of  the 


domain,  and  0  corresponds  to  p   at  the  bottom. 
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The  approximate  steady-state  solutions  to  the  hydrostatic 
primitive  equations  [(2.1)  -  (2.5)]  are: 

V  =  U  =  D  [x(l-e"^cos  C)  +  ye"C  sin  ^]i 

-   y(l-e~^cos  C)  -  xe"^  sin  c)j_,  (2.7) 

a)  =  0, 

<J>  =  $  =  -D2(x2+y2)/2  -  £Dxy  +'  F(p), 
0  =  G(p), 
q  =  H(p), 
where 

f   % 

C  =  (  jy~   )  2 (P  "  P)  >    ^m  ^s  evaluated  at  p  ,  D  and  £  are 

M 
constants  and  where  F(p)  is  such  that  the  hydrostatic  equa- 
tion (2.5)  is  satisfied.   The  functions  G(p)  and  H(p)  are 
arbitrary  functions  of  p  since  w  =  0 .   These  solutions  sat- 
isfy the  boundary  conditions  (2.6)  and  allow  for  inclusion 
of  the  Ekman  boundary  layer.   The  relations  (2,7)  satisfy 
the  steady-state  equations  if  some  small  advection  terms  in 
the  boundary  layer  are  neglected  as  they  usually  are  in 
Ekman  theory. 

The  equations  and  the  model  become  two-dimensional  as  in 
W74  and  Cornelius  [1974]  after  a  subdivision  of  dependent 
variables  as  follows: 

V  =  U(x,y,p)  +  u(y,p,t)i_  +  v(y,p,t)j_ 
u  =  w(y,p,t) 

0  -  e(y,P,t)  (2>8) 

<J>  =  $(x,y)  +  ir(y,p,t) 

q  =  q(y,p,t) 
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It  is  noted  that  all  departures  from  the  approximate 
steady  deformation  solutions  (2.7)  are  assumed  to  be  inde- 
pendent of  x.   If  we  substitute  the  relations  (2.8)  in  the 
component  of  (2.1)  we  obtain: 

|H  +  |_  (uv)  +  |p  (iou)  +  uD(l-eCcos  ?)  -  D[y(l-e~Ccos  C) 

-  xe~C  sin  C]  |~  +  vDe"C  sin  ?  +  coD  |^-  [xe~^(l-cos  c) 

+  ye"^  sin  C]=  fv  +  AM  g»  +  ^  (KM  |»  )  .      (2.9) 

In  this  study  we  will  neglect  the  terms  in  (2.9)  where 
x  appears  and  apply  this  and  the  other  equations  at  x  =  0. 
If  the  initial  u,  v,  w,  9,  r\    and  q  fields  are  independent  of 
x,  the  error  from  this  approximation  will  grow  slowly  and 
will  be  confined  to  the  boundary  layer  since  the  x  dependent 
terms  of  (2.9)  are  zero  outside  the  boundary  layer.   Since 
f rontogenesis  occurs  very  rapidly,  it  is  expected  that  the 
development  of  x-variations  in  the  dependent  variables  would 
have  only  a  small  effect  on  the  resulting  quasi-steady  front. 

In  any  case  these  effects  could  not  be  observed  in  the 
atmosphere,  because  variations  in  the  basic  fields  along 
atmospheric  fronts  would  be  much  more  important. 

When  equation  (2.9)  is  evaluated  at  x  =  0  it  becomes 

3u  3  c  >.  ,  3  ,  ,.  n  3V  PU  3u  3V  ,  -C  „■» 
^r  +  *—  (uv)  +  t^k  (tou)  -  ru  v-  +  TV  7T  ■  v  v  c  ^  sin  C) 
at    3y  v  J         3  P  *■  '  3y       3y      3y  v 

-  V  Ip  (e'?  sin  O  =  fv  ♦  AM  |^  ♦  fg  (*„  f§  ).(2.10) 

where  V  =  -Dy  and  r  =  1-e^  cos  £.   The  j_  component  of  (2.1) 
applied  at  x  =  0  takes  the  form 
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5T  +  §7  ^  +  k   (uv)  +  r  k   (vV)  '  u  ty  e"?  sin  ? 

-  (UV)  |p  f"5  cos  C)  =  -  |f  -  fu  ♦  AM  0   ♦  ,| 

(KMg)   .  C2.ll) 

Similarly  the  following  equation  for  the  departure  of  poten- 
tial temperature  can  be  derived  from  (2.2): 

90,8  r    n-\         9   /-r^    ru  39    .   99,9  rrr      90% 
3t  +  3y  (V8)  +  3?  (U9)  +  rV  37  =  AB  37F  +  3p  (K9  3^  > 

+  QA  +  ^m  +  Hc  t2-12^ 

and  the  following  equation  for  q  from  (2.3) 

|5L  +  |-  (vq)  +  |~  (u)q)  +  TV  |a  =  A   j£f  +  £-  (KM  ^  ) 

9t    9y  v  n^    9p  ^  n^       9y     q  9y2    9p   M  9p  J 

+  MA  -  Mc.         (2.13) 
Equations  (2.4)  and  (2.5)  become 


and 


iZ.  +  1_  =  n  (2  14) 

9y    9P    u'  L^.-Ltj 


15  ■  -L(p)6.  (2.15) 


9p 

The  boundary  conditions  at  the  bottom  and  top  of  the  domain 
respectively,  are 


u  =  v  =  w  =  Ke  ||  =  Kq  fa  =  0    p  =  pQ 

KM  IrT  =  KM  V-    ■  W  =  KA  ¥-    '     K   f^  =  0     P  "  Ptnn 

M  9P     M  9p         0  9p     q  9p       r    rtop 
If  we  define  the  vertical  average  of  a  quantity  as 
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(2.16) 


i  rp° 


P, 
)(dp),  (2.17) 

integrate  the  hydrostatic  equation  (2.15)  with  respect  to  p 
and  remove  the  vertical  mean,  it  follows  that 


P  P0        (2.18) 


7T  -  < 


tt  >  =  [  f       L(p)6dp  -  <    J    L(p)0dp  >]. 


Ptop  Ptop 


Taking  the  vertical  average  of  (2.14)  and  using  the  boundary 
(2.16)  gives 

|y  <  v  >  =  0.  (2.19) 

This  equation  states  that  the  total  mass  flux  in  the  y-direc- 
tion  is  independent  of  y.   In  W72  a  symmetry  argument  was 
used  to  show  that  this  flux  must  be  zero.   This  argument 
does  not  hold  strictly  in  this  experiment,  but  it  can  be 
expected  that  at  large  distances  from  the  frontal  zone,  the 
disturbance  mass  flux  will  vanish.   Thus  it  is  assumed  that 

<  v  >  =  0.  (2.20) 

Equation  (2.11)  can  be  written 

gv  ♦  |_  [vv  -  <vv>]  ♦  fjCuv]  +  |y  (rw-<rv>v) 

-  tt—  [ue~^  sin  c  -  <ue~?  sin  z,>    -    V[w  ■»—  (e~^  cos  c) 
9y  L         *  8p 

-  <w  f-  (e"C  cos  c)>]  ■  fy  (ir-<ir>D  "  f(u-<u>)  +  AM  0- 

+  T-  CKm  I"  )  -  KM  (  1^  )„.„  (2.21) 

3p   M  8p     —   8p   P~P0- 

The  equations  (2.10),  (2.12),  (2.13),  (2.14),  (2.18)  and 
(2.20)  form  a  complete  set  which  can  be  solved  by  a  pure 
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marching  process.   The  finite  difference  equations  conserve 
mean  squares  in  the  advection  terms  and  are  described  by 
Williams  [1967]  .   The  Matsuno  scheme  is  used  everywhere 
except  in  the  arrangement  of  variables  and  the  finite  dif- 
ference approximations  are  the  same  as  those  used  by  Williams 
[1967]  .   In  order  to  close  the  problem,  computational  boun- 
daries must  be  introduced  in  y.   Since  the  disturbance 
velocities  should  die  out  at  a  sufficient  distance  from  the 
axis  of  dilation,  then 

v(±y,p,t)  =  0.  (2.22) 

However,  there  is  appreciable  inflow  across  these  computa- 
tional boundaries  since  V(±Y)  =  ±  DY.  The  quantities  u,  9 
and  q  which  are  advected  across  the  boundaries  must  be  spe- 
cified independent  of  the  interior  values  if  computational 
stability  is  to  be  maintained  (Platzman,  1954).  Thus,  the 
following  boundary  conditions  are  used: 

U[±(Y  +  Ay/2),  p,t]  =  u[±(Y  +  Ay/2),p,0] 

6[±(Y  +  Ay/2),  p,t]  =  9[±(Y  +  Ay/2),p,0] 

8[±(Y  -  Ay/2),  p,t]  =  6[±(Y  -  Ay/2),p,0] 

q[±(Y  +  Ay/2),  p,t]  =  q[±(Y  +  Ay/2),p,0]. 

The  computational  boundaries  y  =  ±Y  are  placed  between  grid 
points  so  that  the  above  conditions  are  actually  applied  at 
y  =  ±(Y  +  Ay/2). 

The  boundary  conditions  were  found  to  be  satisfactory 
except,  that  a  solution  separation  developed  near  the  boun- 
daries.  This  was  controlled  by  using  a  Matsuno  time  dif- 
ferencing scheme  for  stability  and  because  of  its  damping 
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(2.23) 


characteristics.   The  vertical  gridpoints  were  spaced  pro- 


portional to  in   p  using  the  relation  p  =  p   -  p   e 


sk 


where  s  is  Ap/H.   Refer  to  Table  II  for  forecast  levels. 
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III.   INITIAL  CONDITIONS 

The  initial  conditions  of  Cornelius  [1974]  were  modified 

by  raising  the  upper  boundary  to  approximately  the  100-mb 

level,  and  by  simulating  a  "tropopause"  centered  at  about 

the  250-mb  level.   The  Cornelius  [1974]  initial  conditions 

are  developed  as  follows.   The  initial  potential  temperature 

field  in  pressure  coordinates  is  given  by 

c       S?T 

e(y,P,o)  =  eI(po)  exP  [  -|  -£   (l  -  (  2-  )K] 

-    a(2/7r)  arctan(sinh   ay)  (3.1) 

where 


fir  fir 

a    = 


H/LCP^c-aej/Bp)     H/CR/p^Cp/p^^c-aej/ap) 

*  1 


(3.2) 


'M 

In  the  model,  a  is  given  a  constant  median  value  by  solving 

(3.2)  at  p  =  p.,  and  using  the  static  stability  90y/8p  at 
this  level. 

The  initial  u  component  of  the  velocity  is  found  using 
a  form  of  the  hydrostatic  equation  (2.26)  and  the  definition 
of  geostrophic  wind 

»--!$.  (3.3) 

A  form  of  the  thermal  wind  equation  is  found  to  be 

3u  _    1  d      8<J)    1  R  r   p_  vK  dQ  ri    a^ 

Ip~"f8y^p"fp(-poJ   3y  '  l5*4J 

We  again  assume  a  median  pressure  level  p.,  where  u(pM)  =  0, 
and  integrate 

22 


M  M 


to  get 


u(y,P,0)  =  %|f  [(  E-)K  -  (*-*)"]  (3.6) 

7      Fo         *0 

where 

K  =  R/c  . 

p 

From  equation  (4.3),  in  Williams  [1972],  it  is  seen  that 

■sr-  =  -  — —  sech  (ay).  (3.7) 

8y      it         7  * 

Therefore,  the  initial  u  field  is  given  by 

u(y,P,0)  =  -  %   2|«  [(  P_  )K  -  (  ?M  )K]sech  (ay).  (3.8) 

po        *o 

The  initial  divergent  wind  is  determined  as  in  Williams 
[1972]  by  using  the  quasi-geostrophic  equation  in  pressure 
coordinates.   The  stream  function  is  defined  as 

v  =  8ip/8p 
and  (3.9) 

to  =  -  8ip/8y 

and  equation  (4.10)  in  Williams  [1972]  takes  the  form 


d2i>  96 

8pT  +  L^M)("  8p~  ;pM  9y~   fT  ^fMJ  ay 


r  +  L(pMH-  —  )n  7T-  =  TT  l(pm)  Tw  >  C3.10) 

M 


recalling  that 


,  PM  .  K 


L(P"3  =  %   (  P7  ] 

The  equation  (3.10)  is  solved  for  ip  from  the  initial  temper- 
ature field  (3.1)  with  the  technique  of   Ogura  and  Charney 
[1962]  and  ip    is  then  determined  from  (3.9). 
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These  initial  conditions  do  not  satisfy  the  boundary 
conditions  on  u,  v  and  6  at  the  upper  and  lower  boundaries. 
Hence  a  period  of  adjustment  will  be  required  to  form  the 
surface  friction  layer  and  some  oscillations  may  be  observed 
later.   The  fields  of  00  and  tt  are  obtained  through  the  con- 
tinuity and  hydrostatic  equations  (2.14)  and  (2.15). 

The  initial  field  of  specific  humidity  is  calculated 
by  first  determining  saturation  vapor  pressure  e   at  each 
gridpoint  through  use  of  an  integrated  form  of  the  Clausius- 
Clapeyron  equation.   According  to  Hess  [1959] 

Lc    1     1 
es  =  .611  exp  (  g£  (  f-  -  i  ))  (3.11) 

v    o 

where  L   is  the  latent  heat  of  condensation,  R   is  the  gas 
c  '   v        6 

constant  for  water  vapor,  and  T   is  a  constant  of  integra- 
tion.  With  knowledge  of  the  saturation  vapor  pressure, 
the  saturation  specific  humidity  is  obtained  from  the  ap- 
proximate relation  of  Haltiner  and  Martin  [1957] 

.622  e 

q^  = .  (3.12) 

s     p 

Knowledge  of  potential  temperature  at  t  =  0  from  equation 
(3.1),  makes  it  possible  to  calculate  q   at  each  point  using 
equations  (3.11)  and  (3.12).   Defining  q  to  be  a  given 
fraction  of  q   allows  the  testing  of  atmospheres  of  different 
constant  initial  relative  humidity.   The  initial  conditions 
of  Cornelius  [1974]  were  modified  as  follows.   A  tropopause 
was  simulated  centered  at  about  the  250-rnb  level.   This  was 
accomplished  using  a  reference  temperature  at  the  midpoint 
of  the  horizontal  scale  and  calculating  temperatures  such 

24 


that  a  tropopause  with,  its  expected  variation  in  height  was 
determined  and  a  zero  vertical  temperature  gradient  was 
achieved  in  the  lower  stratosphere.   The  tropopause  height 
varied  from  325  to  110  mb  over  the  horizontal  domain  for 
6-j-  (p  )  =  290  K.   The  tropopause  was  simulated  as  follows. 
A  reference  temperature  was  computed  using  equation  (3.1) 
as 

?  P   ~P    -F 

T   ,  =  215  +  a  (  -  )  arctan(sinh  ay)  (  — — —  ) K    (3.13) 
ret  v  tt  J  ' J     K        v>  J      v    J 

*o 

When  the  temperature  in  each  column  reached  this  reference 
temperature,  it  was  held  constant  at  that  value,  as  expected 
in  the  lower  stratosphere.   Temperature  was  then  converted 
back  to  potential  temperature. 
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IV.   CQNVECTIVE  ADJUSTMENT,  LATENT  HEAT  AND  PRECIPITATION 

Two  additional  heating  terms  (Q   and  H  )  were  included 

in  the  moist  thermodynamic  equation  (2.2)  necessitating 

further  amendment  to  W74.   Atmospheric  heating  by  the  dry 

convective  adjustment,  Q  ,  was  included  in  W74  to  eliminate 

a 

any  intense  grid  scale  convection  which  developed  due  to  the 
formation  of  any  dry  unstable  lapse  rates.   This  redistribu- 
tion of  heating  was  accomplished  using  averaging  techniques 
to  ensure  that  the  dry  adiabatic  lapse  rate  was  not  exceeded 
(86/8p  <  0)  . 

Lapse  rates  were  tested  in  each  vertical  column  of  the 
model  as  in  Cornelius  [1974].   Upon  encountering  an  unstable 
layer  (36/8p  >  0) ,  a  cumulative  average  was  taken  and  applied 
to  the  gridpoints  within  this  layer  until  neutral  stability 
(80/8p  =  0)  was  attained.   The  model  was  checked  simultaneously 
for  instabilities  developing  below  the  convectively  adjusted 
layer  and  potential  temperatures  at  these  levels  were  also 
cumulatively  averaged  and  reassigned  the  new  values.   The 
final  result  of  the  technique  is  that  the  heating  is  redis- 
tributed to  simulate  physical  processes  of  the  atmosphere 
and  mathematically  that  99/8p<  0  throughout  each  column  of 
the  model.   A  completely  analogous  scheme  was  used  for  satu- 
rated parcels  in  the  moist  experiment,  only  using  equivalent 
potential  temperatures,  0~E,  in  place  of  potential  temperature 
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The  addition  of  moisture  required  that  large  scale  con- 
densation and  the  subsequent  release  of  latent  heat  be  given 
important  consideration.   Following  Haltiner  [1971],  changes 
in  moisture  and  temperature  due  to  large  scale  condensation 
were  determined  as  follows. 

Where  relative  humidity  exceeds  1001,  the  adjustments  in 
temperature  6T  and  specific  humidity  Sq  are  determined  by 
the  following  equations: 


and 


q  +  6q  =  qs(T  +  6T,p)  (4.1) 


c  6T  =  -  L  6q.  (4.2) 

p        c  n  v    ' 


These  equations  require  that  the  final  state  of  the  air  be 
exactly  saturated  and  that  excess  moisture  be  condensed  iso 
barically,  releasing  the  latent  heat  to  the  air. 

Equations  (4.1)  and  (4.2)  can  be  expressed  in  iterative 

form  as 

8q 
Aq  =  qs(l  -  r)  +  (  -gf  )p  AT,   q  -  rqs         (4.3) 


and 


CAT  =  -LAq.  (4.4) 


Eliminating  AT  from  (4.3)  and  (4.4)  gives 

q  (1  -  r) 

Ac*  =  \ 3qi—  '  (4'5) 

1  +  ^  -7T> 

At  constant  pressure  the  Clausius -Clapeyron  equation  can  be 

expressed  as 

8q        L    q 
r    v 
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Substituting  (4.6)  into  (4.5)  gives 

qs(i  "  r) 
Aq  ■  — th: —  •  .  (4.7) 


1  + 


L2Cls 

c  R  T2 
P  v 


Aq  may  be  now  calculated  from  (4.7)  as  can  the  corresponding 
value  of  AT  from  (4.4).  The  first  estimates  of  the  adjusted 
temperature  and  specific  humidity  are, 

T1  =  T  +  AT    and    q1  =  q  +  Aq.  (4.8) 

These  values  in  turn  are  used  to  obtain  improved  approxima- 
tions.  Iterations  are  continued  until  the  values  of  T'  and 
qf  at  each  level  match  with  sufficient  accuracy  100%  rela- 
tive humidity.   The  changes  are  then  applied  to  the  potential 
temperature  and  specific  humidity  fields  of  the  model  ac- 
counting for  the  term  H   of  equation  (2.2)  and  M  of  equation 
(2.3). 

The  model  is  then  investigated  for  unstable  moist  convec- 
tion and  its  associated  heat  of  condensation.   A  convective 
adjustment  scheme  similar  to  the  dry  one  is  used  on  the  satu- 
rated parcels  so  that  moist  static  stability  is  always  main- 
tained.  A  field  of  equivalent  potential  temperature  is 
formulated  according  to  Holton  [1972]  as 

L  q 

0SE  =  0  exp  (  ~-y  )  (4.9) 

and  the  moist  convective  adjustment  is  designed  to  ensure 
that  the  moist  adiabatic  lapse  rate  is  not  exceeded 
(80SE/3p   <  0). 
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The  heat  release,  Q  ,  as  a  result  of  moist  convective 

^m 

adjustment  is  included  through  knowledge  of  the  change  of 
equivalent  potential  temperature  during  the  adjustment. 

Taking  the  natural  logarithm  of  both  sides  of  equation 

(4.9)  gives 

L  q 
£n  0gE  =  £n  0  ■S-jfc   .  (4.10) 

A  constant  pressure  surface  may  be  assumed  since  the  model 
uses  pressure  coordinates.   Differentiating  both  sides  of 

(4.10)  and  recalling  the  definition  of  potential  temperature 
gives 


dSSE    dT  .  d%    Lc    Lc  %   dT 

eSE  "  t    c  t      t^—  •         t4-1^ 


Solving  equation  (4.6)  for   dq  ,  substituting  into  equation 

(4.11)  and  solving  for  dT  leads  to 

d0c„ 
SE 

dT  =  J q   Lc2 Z-% (4-12) 

q    -\(1    +    — -  - 1 

USE  Tu     R  c^T2     C   T  J 
v  p        p 

where   dO^p  is  the  change  in  0cp  needed  to  adjust  (80ep/3p  = 
0). 

Since  the  atmosphere  remains  saturated  throughout  the 
process  of  moist  convective  adjustment,  the  Clausius-Clapeyron 
equation  is  again  used  to  obtain  the  corresponding  dq  or 
actually  dq  .   Changes  are  again  applied  to  the  potential 
temperature  and  specific  humidity  fields  of  the  model  account- 
ing for  the  terms  QM  of  equation  (2.2)  and  M   of  equation 
(2.3). 
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A  cumulative  amount  of  precipitation  is  then  found  using 
the  relation 

1  r° 

P  =  -  |  J        6qdp  (4.13) 

o 

in  agreement  with  the  changes  in  specific  humidity  due  to 

condensation  and  moist  convective  adjustment.   Precipitation 

is  expressed  in  inches  of  rainfall  in  the  model. 
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V.   NUMERICAL  SOLUTIONS 

All  numerical  results  to  be  shown  use  the  following 
values  for  the  constants: 

At  =  540  sec         p   =  100  cb  pw  =  55  cb 

*o  rM 

Ap/H  =  .01978  Y  =  1800  KM  Ay  =  60  KM 

g  =  9.81  m  sec"2      f  =  10  "*  sec"1       D  =  10"1  sec"1 

0   =  300  K       6T(p  )  =  280,290,300  K  c   =  1003  joules 
0                *   °  P       Kg"1  Kg1 

R  =  461  joules       R  =  287  joules  Kg"1  K*1 


v 


Kg   K 


k  =  R/c  L   =  2.5xl06  joules 

K,"1 


a  =  12.56  K         RH  =  70% 

90  80 

-3^-  =  5  K  Km"1      (-^i.)    =  -.53448  K  cb" 
3z  \  3p  P, 


M 
"m    "o  "  'q 

i 

'q 


A„  =  A„  =  3xl04m2sec"  A_  =  100000  m2sec" 

C   =  0,  cb2sec_1       CmB  =  0,  .001409,  .0007  cb2sec 


(which  corresponds  to  0,  10,  5  m2sec    in  height  coordinates) 

CQB  =  0,  .0007,  .0001761  cb  sec    (which  corresponds  to 

0,  5,  1.25  m2sec    in  height  coordinates) 

where  more  than  one  value  means  that  the  constant  was  varied 

in  different  experiments.   Whenever  possible,  the  values  are 

the  same  as  those  used  in  W74. 

The  values  of  vertical  diffusion  of  moisture  C   and  heat 

m 

Cfi  were  varied  in  the  experiments  to  compensate  for  density 
variation  in  the  atmosphere.   This  is  equivalent  to  a  constant 
coefficient  in  the  z  coordinate  system  and  is  accomplished 
using  the  following  formulae: 
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and 

Ce  -  CeB(l  -  P/P0)2  (5.2) 

in  this  section  numerical  f rontogenesis  solutions  with  mois- 
ture, friction  and  a  stratosphere  are  examined  at  t  =  4  days 
and  compared  to  the  corresponding  dry  atmosphere. 

As  can  be  seen  from  Table  I,  11  experiments  were  performed 
Experiments  1  and  2  contain  only  horizontal  turbulent  diffu- 
sions of  momentum,  heat  and  moisture  while  Experiments  3 
through  11  contain  both  horizontal  and  vertical  turbulent 
diffusions  of  these  quantities. 

In  Table  I,  9T(p  )  is  the  reference  potential  temperature 
(K) .   The  condensation  interval  is  5  time  steps  between  suc- 
cessive applications  of  condensation  and  moist  convective 
adjustment  to  the  model  in  the  moist  experiments.   Horizon- 
tal turbulent  diffusion  coefficients  and  initial  boundary 

values  of  vertical  turbulent  diffusion  coefficients  (C  -r, , 

*■  mB' 

CfiR,  C  )  are  relative  humidity  (percent)  and  initial  static 
stability  (K/cb) . 

Experiments  were  performed  with  a  reference  potential 
temperature  (G\.  (p  )),  of  280K,  290K  and  300K  to  simulate 
physically  realistic  situations.   Changes  to  the  dry  atmo- 
sphere with  different  reference  temperatures  were  almost 
negligible,  but  due  to  the  exponential  function  in  the 
Clausius.-Clapeyron  equation  important  differences  were  ob- 
served in  the  moist  atmosphere. 
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Experiments  were  also  conducted  using  different  boundary- 
values  of  vertical  turbulent  diffusion  coefficients  of  momentum 
and  heat. 

The  final  experiment  included  a  turbulent  mixing  coeffi- 
cient which  reduced  the  values  of  CmP  and  C^  linearly  to  .1 

inn      QB 

of  their  original  values  between  approximately  the  93-  and 
84-cb  levels  in  five  equal  steps. 

The  adjusted  values  of  vertical  turbulent  diffusion  co- 
efficients of  momentum  C   and  heat  CA  were  then  varied  as 

m  6 

described  by  equations  (5.1)  and  (5.2). 

Primary  emphasis  on  analyses  of  experimental  results  was 
on  contour  plots  from  the  IBM  360/67  computer.   These  results 
obtained  for  all  11  experiments  are  analyses  at  t  =  4  days 
of  the  fields  of  u,  q,  0  and  w.   The  entire  vertical  plane 
is  shown  from  100  to  10  cb  in  these  plots.   However,  the 
outer  portion  of  the  domain  between  |y|  =  1200  km  and  |y|  = 
1800  km  is  not  shown. 

Figures  1-3  show  the  initial  fields  of  u,  8  and  w  at 
reference  potential  temperature  of  290  K  used  in  the  follow- 
ing experiments.   Contour  plots  for  experiments  1  and  2  are 
shown  as  figures  4-10.   Figures  4  and  5  show  distributions 
of  the  rotational  wind,  u  for  corresponding  dry  and  moist 
experiments,  respectively.   Features  of  u  at  all  levels  dif- 
fer quite  significantly  from  dry  to  moist  experiments.   There 
is  a  marked  increase  in  horizontal  and  vertical  shear  in  the 
easterly  flow  in  the  moist  experiment  in  the  vicinity  of  the 
front  near  the  surface.   Horizontal  and  vertical  wind  shear 
is  also  more  intense  at  low  middle  to  middle  levels  in  the 
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moist  experiment.   It  is  noteworthy  that  this  property  is 
dominant  throughout  the  moist  solutions.   In  the  moist  ex- 
periment the  jet  core  increases  in  magnitude  while  shear 
decreases  to  zero  at  a  lower  level  in  the  northern  region 
near  the  tropopause.   Figures  6  and  7  contrast  values  of 
e  throughout  the  domain.   Noted  is  the  warmer  air  in  middle 
levels  in  the  frontal  zone.   A  slightly  more  intense  gra- 
dient is  seen  near  the  surface  with  little  change  in  tilt 
of  the  front.   All  this  is  attributed  to  the  latent  heat 
release  in  the  moist  experiment  2. 

Vertical  motions  for  the  experiments  are  investigated  in 
Figures  8  and  9.   Upward  vertical  motions  occupy  a  much 
larger  zone,  extend  much  higher  and  occur  closer  to  the 
frontal  zone  in  the  warm  air  mass  of  the  moist  experiment. 
Significant  downward  motion  occurs  between  |Y|  =  1200  and 
|Y|  =  1800  km  (not  shown  in  the  figures). 

Figure  10  shows  the  field  of  specific  humidity  for  ex- 
periment number  2.   Maximum  precipitation  is  .37  inches. 
Figures  11  through  17  compare  the  dry  and  moist  fields  of 
u,  0  and  w  and  show  the  q  fields  of  experiments  3  and  4. 
Figures  11  and  12  show  a  small  increase  in  the  zone  of 
easterly  component  in  the  lower  levels  of  the  moist  experi- 
ment.  Vertical  and  horizontal  shear  are  also  seen  to  in- 
crease as  in  experiments  1  and  2.   Also  noted  is  a  much 
larger  jet  core  in  the  moist  experiment,  with  a  slight 
northward  shift.   Noteworthy  is  a  significant  decrease  in 
magnitude  of  the  easterly  component  of  u  near  the  surface 
from  that  of  experiments  1  and  2.   Figures  13  and  14  exhibit 
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the  same  contrast  in  potential  temperature  as  that  of  ex- 
periments 1  and  2.   It  is  noted  that  the  air  at  low  levels 
and  near  the  tropopause  is  colder  in  the  moist  experiment 
with  vertical  turbulent  diffusion  added  (exp.  4).   This  is 
also  true  in  the  dry  experiments.   This  is  attributed  to 
the  increased  amount  of  heat  loss  due  to  the  vertical  tur- 
bulent diffusion  of  momentum  and  heat  offsetting  the  release 
of  latent  heat  in  the  moist  experiments.   The  air  at  mid 
levels  in  the  frontal  zone  is  warmer  in  the  moist  experi- 
ment as  seen  before.   Figures  15  and  16  show  a  marked  change 
in  vertical  motion  Jbetween  dry  and  moist  experiments.   A 
large  zone  of  upward  motion  is  seen  near  the  frontal  zone 

-  3 

with  maximum  values  exceeding  3x10    mb/sec  near  the  85-cb 
level.   Figure  17  shows  the  specific  humidity  of  experiment 
4.   Noteworthy  is  the  dramatic  shift  to  the  south  in  specific 
humidity  following  the  frontal  position. 

Figures  19  through  38  compare  the  moist  and  dry  fields 
of  u,  0  and  w  and  show  the  specific  humidity  of  the  moist 
experiments  using  boundary  values  of  vertical  turbulent  dif- 
fusion of  heat  and  momentum  of  .0007  cb2sec    (5  m2sec 
in  height  coordinates)  at  reference  potential  temperatures 
of  290,  280  and  300  K. 

Figures  18  and  19  contrast  the  u  field  of  experiments  5 
and  6  at  "6 T  (p  )  =  290  K.   Noteworthy  is  the  increase  in  zone 
of  maximum  easterly  component  of  the  wind  at  low  levels  and 
the  more  intense  horizontal  and  vertical  shear  in  low  levels 
over  experiments  3  and  4.   A  slight  northward  shift  and  in- 
crease in  size  of  the  jet  core  is  noted  in  the  moist  experiment 
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Figures  20  and  21  show  little  difference  in  potential  tem- 
perature fields,  although  the  air  in  the  frontal  zone  in 
middle  levels  is  warmer  in  the  moist  experiment  due  latent 
heat  release.   There  is  a  marked  increase  in  tilt  of  the 
front  over  previous  experiments,  due  to  more  normal  north- 
ward position  of  the  front  at  the  surface.   Figures  22  and 
23  contrast  the  vertical  motion  of  experiments  5  and  6  show- 
ing a  maximum  in  upward  motion  near  the  850-mb  level  con- 
sistent with  previous  experiments.   Figure  24  shows  the 
specific  humidity  field  which  differs  slightly  more  in  the 
same  manner  as  experiment  4  differed  from  2.   Maximum  pre- 
cipitation is  .4  in. 

Figures  25  and  26  show  more  intense  horizontal  and  ver- 
tical shear  apparent  in  the  easterly  component  of  the  moist 
experiment  at  9" T  (p  )  =  280  K.   The  jet  core  is  reduced  in 
size  and  magnitude  from  experiments  5  and  6  and  the  hori- 
zontal shear  decreases  significantly  near  the  tropopause. 
It  should  be  noted  here  that  the  tropopause  in  experiments 
7  and  8  is  lower  than  the  tropopause  in  previous  and  subse- 
quent experiments  due  to  the  lower  reference  potential  tem- 
perature used.   Figures  27  and  28  show  little  difference 
between  moist  and  dry  fields  of  potential  temperature  except 
for  the  warmer  air  in  mid  levels  in  the  frontal  zone  of  the 
moist  experiment  due  latent  heat  release.   Noteworthy  in 
these  experiments  7  and  8  is  the  zone  of  low  stability  and 
apparent  break  in  the  tropopause  in  the  northern  part  of  the 
region.   Figures  29  and  30  show  a  similar  contrast  between 
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moist  and  dry  fields  of  vertical  motion  as,  seen  in  previous 
experiments  except  that  the  contrast  is  not  as  dramatic  as 
seen  before.   Figure  31  shows  the  q  field  is  reduced  con- 
siderably over  previous  experiments  due  to  the  lower  refer- 
ence potential  temperature  used  (280  K) .   Maximum  precipita- 
tion reduced  to  .15  inches. 

Figures  32  through  40  show  a  much  higher  tropopause 
due  to  the  larger  reference  potential  temperature  (300  K) . 
A  dramatic  change  in  the  u  field  is  seen  in  Figures  32  and 
33  between  dry  and  moist  fields.   Especially  noteworthy 
is  the  very  strong  jet  core  and  much  more  intense  horizon- 
tal and  vertical  wind  shear.   Also  note  a  large  shift  in 
the  jet  to  the  north  which  was  not  quite  so  obvious  in  pre- 
vious experiments.   Figures  34  and  35  also  show  a  marked 
contrast  between  dry  and  moist  potential  temperature  fields. 
The  air  is  much  warmer  in  the  upper  levels  near  the  tropo- 
pause due  to  the  latent  heat  release.   Thermal  gradients 
along  the  frontal  zone  are  also  stronger  in  the  mid  and  up- 
per levels.   Strong  heating  in  upper  levels  is  apparent. 

-  3 

Vertical  motions  (Figures  36  and  37)  exceed.  3x10   mb/sec 
in  the  warm  moist  air.   The  contrast  is  similar  to  that  seen 
in  previous  experiments.   The  moisture  field  (Figure  38)  is 
greatly  intensified  over  previous  experiments  and  specific 
humidities  are  in  excess  of  19  g/kg  at  low  levels  to  the 
south.   This  is  no  doubt  due  to  the  larger  quantities  of 
latent  heat  release  observed  in  this  experiment,  and  the 
larger  initial  reference  potential  temperature.   Maximum 
precipitation  increased  to  .71  inches. 
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Figures  39  and  42  show  experiment  11  which  included  the 
turbulent  mixing  coefficient  into  experiment  4.   The  only 
differences  noted  were  a  small  pocket  of  lower  specific 
humidity  (Figure  42)  at  the  top  of  the  mixing  level  and  a 
decreased  zone  of  maximum  upward  motion  in  the  warm  moist 
air  of  the  mixed  experiment  (Figure  41) .   From  the  above 
we  conclude  that  vertical  turbulent  mixing  has  no  effect 
on  mid  and  upper  level  deformation  fields. 

Comparison  of  the  results  of  the  previous  experiments 
with  similar  experiments  of  Cornelius  [1974]  revealed  the 
following : 

1.  The  westerly  wind  field  extended  much  lower  while 
the  easterly  wind  field  was  compressed  below  the  85  cb 
level.   The  zero  wind  level  was  also  observed  to  be  lower 
in  the  expanded  atmosphere  with  tropopause  included.   This 
is  true  in  both  moist  and  dry  experiments.   The  jet  was  also 
lower  and  more  apparent. 

2.  The  pockets  of  warm  air  were  apparent  in  the  frontal 
zone  as  in  Cornelius  [1974]  ,  but  frontal  tilt  was  somewhat 
less  after  4  days  in  the  experiments  conducted  in  this  study. 

3.  Significant  downward  vertical  motion  occurred  outside 
the  zone  observed  in  both  studies  although  larger  values  are 
evident  in  the  studies  of  Cornelius  [1974] . 

4.  Specific  humidity  fields  were  similar  in  both  studies 


38 


VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  numerical  model  of  Williams  [1974]  is  modified  in 
this  study  to  include  moisture,  a  tropopause  and  lower  stra- 
tosphere and  variable  vertical  turbulent  diffusion  of  heat 
and  momentum.   This  leads  to  a  more  realistic  situation  and 
consequently  more  realistic  solutions.   Intensification  of 
the  baroclinic  zone  is  seen  to  occur  in  upper  levels,  al- 
though no  effect  is  apparent  on  the  surface  due  to  the  in- 
clusion of  moisture.   The  effects  of  vertical  turbulent 
diffusion  is  apparent  throughout  the  frontal  zone  through 
increased  tilt,  reduced  heating  and  reduced  values  of  wind- 
speed.   The  expanded  boundary  and  inclusion  of  the  tropo- 
pause along  with  a  conversion  to  In   p  vertical  spacing 
tended  to  shift  all  features  to  a  lower  level  in  the  domain. 
Jet  cores  were  generally  larger  and  vertical  motions  in- 
creased in  the  warm  moist  air  near  the  frontal  zone  due  to 
the  heating.   The  moist  models  illustrate  an  exponential 
thermal  dependence  in  the  various  experiments  at  the  reference 
potential  temperatures  of  280,  290  and  300  K. 

In  future  studies  using  this  model,  attention  should  be 
directed  toward  the  use  of  still  more  realistic  boundaries, 
such  as  a  higher  upper  boundary.   A  translational  motion 
should  be  added  to  the  system  so  that  warm  and  cold  frontal 
effects  could  be  investigated.   Finally,  an  improved  convec- 
tive  scheme  such  as  that  of  Arakawa  and  Schubert  [1973] 
could   be  applied  to  the  model. 
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Table  II.   Forecast  Levels  (cb) 


LEVEL 

p(cb) 

p(cb) 

p(cb) 

1 

95.86 

19 

44.81 

37 

20.95 

2 

91.9 

20 

42.96 

38 

20.08 

3 

88.1 

21 

41.18 

39 

19.25 

4 

84.45 

22 

39.48 

40 

18.25 

5 

80.96 

23 

37.85 

41 

17.69 

6 

77.61 

24 

36.2  8 

42 

16.96 

7 

74.4 

25 

34.78 

43 

16.26 

8 

71.32 

26 

33.34 

44 

15.59 

9 

68.37 

27 

31.96 

45 

14.94 

10 

65.54 

28 

30.64 

46 

14.32 

11 

62.83 

29 

29.37 

47 

13.73 

12 

60.23 

30 

28.16 

48 

13.16 

13 

57.54 

31 

26.99 

49 

12.62 

14 

55.35 

32 

25.88 

50 

12.10 

15 

53.06 

33 

24.81 

51 

11.60 

16 

50.87 

34 

23.78 

52 

11.12 

17 

48.76 

35 

22.80 

53 

10.66 

18 

46.75 

36 

21.85 

54 

10.22 
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